Frustrated couplings between alternating spin-i chains in AgV0As04 



O 

< 

0\ 



Alexander A. Tsirlin/'B Ramesh Nath/'^ Jorg Sichelschmidt,^ 
Yurii Skourski,'^ Christoph Geibel/ and Helge Rosner^'0 

^Max Planck Institute for Chemical Physics of Solids, Nothnitzer Str. 40, 01187 Dresden, Germany 
^School of Physics, Indian Institute of Science Education and Research, Trivandrum-695016 Kerala, India 
^Dresden High Magnetic Field Laboratory, Helmholtz-Zentrum Dresden-Rossendorf, 01314 Dresden, Germany 

We report on the crystal structure and magnetic behavior of the spin-i compound AgV0As04. 
Magnetic susceptibility, high-field magnetization, and electron spin resonance measurements identify 
AgVOAs04 as a gapped quantum magnet with a spin gap A ~ 13 K and a saturation field fioHs — 
48.5 T. Extensive band structure calculations establish the microscopic magnetic model of spin chains 
with alternating exchange couplings J ~ 40 K and J' ~ 26 K. However, the precise evaluation of 
the spin gap emphasizes the role of interchain couplings which are frustrated due to the peculiar 
crystal structure of the compound. The unusual spin model and the low energy scale of the exchange 
couplings make AgV0As04 a promising candidate for an experimental investigation of Bose-Einstein 
condensation in high magnetic fields. 

PACS numbers: 75.50.-y, 75.30.Et, 71.20.Ps, 61.66.Fn 
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I. INTRODUCTION 



Quantum spin systems exhibit diverse and unusual 
phenomena in high magnetic fields While the bulk of 
the experimental and theoretical work refers to systems 
of weakly coupled spin dimers, similar physics is expected 
in bond-alternating (dimerized) spin-i chains that like- 
wise feature a singlet ground state with a spin gap. The 
couplings between the dimers or dimerized chains tend 
to reduce the spin gap, and heavily affect the high- 
field behavior. For example, Bose-Einstein condensa- 
tion (BEG) of magnons, which is a general paradigm 
for field-induced antiferromagnetic (AFM) order in dimer 
magnets, exhibits an exotic two-dimensional (2D) regime 
in BaCuSi206, where interdimer couplings are partially 
frustrated;^ The abundance of intriguing theoretical pre- 
dictions, especially for gapped frustrated systems f2r^ mo- 
tivates the search for low-dimensional spin-gap magnets 
with unusual and potentially frustrated geometry of the 
exchange couplings. 

Efficient experimental studies require materials with 
weak exchange couplings that allow exploring the full 
temperature- vs-field phase diagram up to the saturation 
field Hs- Direct connections between the polyhedra of 
magnetic atoms often lead to strong exchange couplings 
on the order of 100 K. Resulting saturation fields well 
exceed 100 T and remain unfeasible for present-day ex- 
perimental facilities. In structures with non-magnetic 
polyanions, such as phosphates, silicates, or borates, 
the superexchange is dominated by M-O-O-M pathways 
(here, M is a transition metal). The long pathways re- 
quire a specific arrangement of atomic orbitals and usu- 
ally lead to weak exchange couplings on low-dimensional 
spin lattices. For example, phosphates are known 
as a playground for studying spin-i frustrated square 
lattices .i^^-iS Further on, symmetry restrictions on mag- 
netic orbitals induce a variety of one-dimensional (ID) 
spin lattices in vanadium compounds with 2D^^'^^ or even 



three-dimensional (3D)i^"— crystal structures. 

In the following, we present an interesting implementa- 
tion of the quasi-lD spin lattice in the 3D crystal struc- 
ture of AgV0As04. An extensive experimental and com- 
putational study identifies this compound as a system of 
alternating spin chains with frustrated interchain cou- 
plings. The saturation field of about 50 T enables ex- 
perimental access to the full temperature-vs-field phase 
diagram, and makes AgV0As04 an appealing material 
for high-field studies that could probe the BEG physics. 
The outline of the paper is as follows. After reporting 
methodological aspects in Sec. |lTl we present the crys- 
tal structure from high-resolution synchrotron powder 
diffraction (Sec. lIIII) . the experimental study by magneti- 
zation and electron spin resonance (ESR) measurements 
(SecHYl), and the microscopic magnetic model from den- 
sity functional theory (DFT) band structure calculations 
(Sec. |V]). The experimental and computational results 
are compared in Sec. IVII that also provides an outlook 
for future experiments. 



II. METHODS 

Powder samples of AgV0As04 were prepared by an- 
nealing the stoichiometric mixture of Ag20, VO2, and 
AS2O5 in an evacuated silica tube at 600 °G for 24 
hours. To avoid the hydration of AS2O5, the reagents 
were handled in an Ar-filled glove box. According to x- 
ray diffraction (XRD) data collected with Huber G670 
Guinier camera (GuK^i radiation, 29 = 3 — 100° an- 
gle range, image- plate detector), the samples contained 
the AgV0As04 phase and a minor impurity of metallic 
silver [the weight fraction of Ag is 1.5(1) %, as found 
from the Rietveld refinement]. The slight reduction of 
silver is likely caused by a ready decomposition of Ag20 
at elevated temperatures, especially in vacuum. How- 
ever, our attempts to prepare AgV0As04 via annealing 
in air failed due to the partial oxidation of V^'*. Since 
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TABLE I. Lattice parameters of AgV0As04 and refinement residuals Ri/Rp. The standard devi- 
ations are taken from the Rietveld refinement. 



r(K) 


a (A) 


MA) 


c(A) 


13 (deg) 


Ri/Rp 


20 


6.70096(1) 


8.84034(2) 


7.26478(2) 


115.180(1) 


0.028/0.090 


300 


6.71576(1) 


8.84882(2) 


7.28481(2) 


115.282(1) 


0.033/0.068 



the metallic silver shows temperature-independent mag- 
netic susceptibility, this impurity does not influence any 
of the results presented below. Note that the additional 
annealing of the AgV0As04 samples above 600 °C led 
to the melting and decomposition of the compound. 

High-resolution XRD data for structure refinement 
were collected at the ID31 beamline of European Syn- 
chrotron Radiation Facility (ESRF) with a constant 
wavelength of about 0.4 A. The signal was measured by 
eight scintillation detectors, each preceded by a Si (111) 
analyzer crystal, in the angle range 26* = 1 — 40°. The 
powder sample was contained in a thin-walled borosili- 
cate glass capillary with an external diameter of 0.5 mm. 
To achieve good statistics and to avoid the effects of the 
preferred orientation, the capillary was spun during the 
experiment. The program JANA2000 was used for the 
structure refinement. -'^^ 

For magnetization measurements in fields below 5 T 
and in the 2 — 380 K temperature range, we used 
the laboratory Quantum Design SQUID magnetome- 
ter (MPMS 5 T). At higher fields, the magnetization 
isotherm was collected at 1.4 K in the Dresden High Mag- 
netic Field Laboratory using a pulsed magnet. Details 
of the experimental setup are described elsewherCfii The 
perfect coincidence of the data collected during the uprise 
and decrease of the field indicated the lack of irreversible 
effects upon the magnetization of the sample. 

To simulate magnetic susceptibility and high-field 
magnetization curves, we used the ALPS simula- 
tion packagei^ and the quantum Monte-Carlo di- 
rected loop algorithm in the stochastic series expansion 
representation.^^ The simulations were performed for fi- 
nite lattices with periodic boundary conditions. The size 
of the finite lattice was fixed to L = 120 for ID systems 
and L X L with L = 24 for 2D systems. The convergence 
with respect to the cluster size and to the aspect ratio 
(in the case of 2D clusters) was carefully checked. 

The ESR experiments were performed with a commer- 
cial CW spectrometer for temperatures between 10 K 
and 294 K by using a "^He flow- type cryostat. As a func- 
tion of an external, static magnetic field B, the resonance 
shows up as absorbed power P of a transversal magnetic 
microwave field at a frequency of 9.4 GHz ("X-band") or 
34 GHz ("Q-band"). To improve the signal-to- noise ra- 
tio, a lock-in technique was used by modulating the static 
field which yields the derivative of the power absorption 
dP/dB. 

Scalar-relativistic DFT band structure calculations 
were performed within the Perdew-Wang^'^ local density 



approximation (LDA) for the exchange-correlation po- 
tential. The pre-optimizcd basis set comprised atomic- 
like local orbitals of the FPLO scheme (FPL08 . 50)^ The 
calculations were based on the low-temperature crystal 
structure refined from the XRD data. The use of the 
room-temperature structural data led to only minor dif- 
ferences in individual exchange couplings. 

For the LDA calculations, the crystallographic unit 
cell and the /c-mesh of 4096 points (1170 in the irre- 
ducible part) were used. The LDA band structure was 
mapped onto a multi-orbital tight-binding (TB) model, 
whereas correlations effects in the V 3d shell were treated 
within the Hubbard model. Exchange couplings were de- 
rived in the limit of strong correlations (model approach) . 
Alternatively, the correlation effects were treated in a 
mean-field approximation by local spin density approxi- 
mation (LSDA)-|-[/ calculations for supercells (applying 
a /c-mesh of 144 points). Total energies for a number 
of coUinear spin configurations were mapped onto the 
Heiscnberg model (supercell approach). Further details 
of the computational procedure are given in Sec. El 



III. CRYSTAL STRUCTURE 

The room-temperature x-ray powder pattern of 
AgV0As04 was indexed in a monoclinic unit cell with 
lattice parameters a ~ 6.712 A, 6 ~ 8.849 A, c ~ 7.285 A, 
and /3 ~ 115.28°. The refiection conditions 0A:0, k = 2n 
and hOl, I = 2n identified the P2i/c space group. The 
close correspondence to NaV0As04 (P2i/c, a ~ 6.65 A, 
b ~ 8.71 A, c ~ 7.22 A, and /3 ~ 115.2°)^ suggests a 
similarity of the crystal structures, as further confirmed 
by the successful Rietveld refinement, ^'^ Atomic displace- 
ment parameters (ADPs) of oxygen atoms were con- 
strained (refined as a single parameter) due to the sub- 
stantial difference in the scattering powers of Ag and O. 

Tables U and [iTl list lattice parameters and refined 
atomic positions for AgV0As04 at 20 and 300 K. Upon 
cooling, the structure shows an almost uniform shrink- 
age along all three directions. The atomic positions 
are nearly unchanged. We also did not observe any 
non-indexed refiections or reflection splittings that could 
indicate a lower symmetry or a superstructure.^* The 
noteworthy feature is the relatively high ADP of Ag at 
300 KM- At 20 K, the ADP of silver is, however, re- 
duced and becomes comparable to ADP's of other atoms. 
Therefore, the high ADP at 300 K should be assigned to 
thermal vibrations and does not indicate an intrinsic dis- 



FIG. 1. (Color online) Left: crystal structure of AgV0As04 with curved lines separating the spin chains. Filled circles show 
the Ag atoms. Middle: a sketch of the spin model; thick (dark-yellow) line denotes the four-member frustrated unit. Right: 
the frustrated unit superimposed on the crystal structure. Open circles depict vanadium positions/sites of the spin lattice. 



TABLE II. Atomic positions and isotropic atomic displace- 
ment parameters {Uiso, in units of 10"'^ A'^) for AgV0As04 
at 20 K (first set of lines) and 300 K (second set of lines). 
All atoms occupy the general 4e position of the P2\/c space 
group. 



Atom 


x/a 


y/b 


z/c 




Ag 


0.2418(1) 


0.0846(1) 


0.7204(1) 


0.13(1) 




0.2424(1) 


0.0847(1) 


0.7225(1) 


1.49(2) 


V 


0.7441(2) 


0.2376(2) 


0.0300(2) 


0.07(2) 




0.7445(2) 


0.2383(1) 


0.0308(2) 


0.39(3) 


As 


0.2494(2) 


0.0606(1) 


0.2466(1) 


0.11(1) 




0.2496(2) 


0.0605(1) 


0.2468(1) 


0.43(2) 


0(1) 


0.8435(8) 


0.0487(5) 


0.6140(6) 


0.00(4) 




0.8518(8) 


0.0493(5) 


0.6216(6) 


0.38(4) 


0(2) 


0.0466(8) 


0.1722(5) 


0.0750(6) 


0.00(4) 




0.0457(8) 


0.1715(5) 


0.0760(6) 


0.38(4) 


0(3) 


0.7488(9) 


0.1912(4) 


0.2498(7) 


0.00(4) 




0.7469(10) 


0.1895(4) 


0.2489(8) 


0.38(4) 


0(4) 


0.6371(8) 


0.0445(5) 


0.8797(6) 


0.00(4) 




0.6354(8) 


0.0439(5) 


0.8762(6) 


0.38(4) 


0(5) 


0.4339(8) 


0.1845(5) 


0.4077(6) 


0.00(4) 




0.4359(8) 


0.1819(5) 


0.4082(6) 


0.38(4) 



order. 

The crystal structure of AgV0As04 (Fig. [T] and Ta- 
ble llllj) is essentially similar to that of NaV0As04. Vana- 
dium atoms form distorted octahedra with the short 
vanadyl bond of about 1.64 A. The V~0 separations 
in the equatorial plane of the octahedron amount to 
1.99-2.01 A, while the longest V-0 distance of 2.14 A is 
found opposite to the vanadyl bond. Arsenic atoms form 
nearly regular ASO4 tetrahedra with As-0 distances of 
1.69— 1.71 A. Vanadium octahedra share 0(3) corners to 
form chains along the c direction. The ASO4 tetrahedra 
link the neighboring chains into a 3D framework. 

Magnetic properties of the NaV0As04-type com- 
pounds have not received a detailed investigation. The 



TABLE III. Interatomic distances (in A) in the AgV0As04 
structure at 20 K. 



Ag 


5-0(l) 


2.463(4) 


V-O(l) 


2.010(4) 






2.536(4) 


V-0(2) 


1.996(5) 


Ai 


!-0(2) 


2.505(4) 


V-0(3) 


1.636(6) 




;-0(3) 


2.446(4) 




2.144(6) 


Ai 


!-0(4) 


2.423(5) 


V-0(4) 


1.988(4) 






2.897(4) 


V-0(5) 


2.003(5) 


Ai 


;-o(5) 


2.487(4) 


As-O(l) 


1.698(5) 








As-0 (2) 


1.712(4) 








As-0 (4) 


1.697(5) 








As-0(5) 


1.693(4) 



magnetic susceptibility of NaV0P04 was tentatively 
analyzed within the uniform-chain model, based on 
the uniform chains of the VOe octahedra ("structural 
chains") running along the c direction. The reference to 
other V+^ compounds, however, does not support this 
interpretation. A distinctive feature of V"'"'' is the forma- 
tion of a short bond to oxygen that, in turn, places the 
unpaired electron of vanadium on the dxy orbital (here, 
z aligns with the short bond and the structural chain). 
This arrangement of the magnetic orbital precludes the 
V-O-V superexchange along the structural chains, and 
makes the relationship between the crystal structure and 
the spin lattice rather complexj ^^i^^i^^ In Sec.|Vl we apply 
extensive DFT calculations in order to find the correct 
microscopic model. Similar to other V^'' compounds, the 
structural chains do not represent the direction of the 
leading exchange in AgV0As04, although the magnetic 
behavior is quasi-lD. 
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O Experiment (0.5 T) 
alt. chain (ID) 




Temperature (K) 

FIG. 2. (Color online) Magnetic susceptibility of AgVOAs04 
measured in the applied field of 0.5 T (circles) and the fit with 
the alternating-chain model [Eq. 0] (solid line). The inset 
shows the Curie- Weiss fit [Eq. ([T|, dashed line]. 



IV. MAGNETIC PROPERTIES 
A. Magnetization 

Temperature dependence of the magnetic susceptibil- 
ity (x) is shown in Fig.[2]and reveals paramagnetic Curie- 
Weiss behavior above 100 K. At low temperatures, we 
find a maximum around 20 K followed by an increase in 
the susceptibility below 6 K. The susceptibility maximum 
is a characteristic feature of a low-dimensional and/or 
frustrated system with predominantly antiferromagnetic 
(AFM) exchange couplings. The low-temperature upturn 
can be ascribed to the paramagnetism of defects and im- 
purities, as further confirmed by the suppression of this 
upturn in magnetic field. ^'^ Note that we tried to vary 
the preparation procedure in order to improve the sam- 
ple quality, but the low-temperature feature remained 
nearly unchanged. 

The data in the paramagnetic region (above 100 K— ) 
are fitted with the Curie- Weiss law (inset of Fig. [2): 



X = Xo + 



c 



T + e' 



(1) 



where xo accounts for temperature-independent Van 
Vleck and diamagnetic contributions, C stands for the 
Curie constant, and 9 is the Curie- Weiss temperature. 
We find Xo = -2.0(1) x lO""* emu/mol, C = 0.363(2) emu 
K/mol, and ^ 20.6(6) K. The C value yields the effec- 
tive magnetic moment of fics = 1.70(1) fiB that closely 
approaches the spin-only value of 1.73 fiB for V"*''^ and 
corresponds to g = 1.96 in good agreement with the ESR 
result of 1.947 (see below). The positive Curie- Weiss 
temperature 9 identifies leading AFM couplings. 

To fit the susceptibility maximum, a low-dimensional 
spin model should be applied. Although not obvious at 
this stage, we use the ID model of an alternating spin 
chain, as found from the microscopic analysis (Sec. \V\ . 
From a phenomenological point of view, this model can 



be justified by the high-field magnetization curve (Fig.|3]) 
that indicates a spin gap in AgV0As04. The gapped 
nature of the excitation spectrum immediately rules out 
the gapless uniform-chain model. Therefore, the naive 
structure-based assignment of the spin model (Sec. IIIII 
and Ref. apparently fails. 

Since the data at low temperatures are dominated by 
the paramagnetic impurity contribution, we fit the sus- 
ceptibility curve with 



X = Xo + 



imp 



T + 9u 



Xid{J,J') 



(2) 



where Xo has the same meaning as in Eq. ([T]), the 
Ch-np/{T + 9imp) term accounts for the low-temperature 
impurity contribution with ^imp being an effective in- 
teraction between the impurity spins, and xiJ^ J') 
the susceptibility of the alternating spin-i chain with 
interactions J and J' (at this point, we do not dis- 
cuss the relation of these couplings to the crystal 
structure))^ Overall, Eq. ([2]) has six variable parame- 
ters: xo,C'iinp,6'imp,57'^, and J'. The fit (Fig. [1]) yields 
XO = -1.9(1) X 10"^ emu/mol, dmp = 0.027(1) emu 
K/mol (7 % of spin-i paramagnetic impurities), 6'iinp = 
1.3(1) K,g = 1.87(1), J = 41.8(1) K, and J' = 25.8(1) K. 
The fitted value is lower than g = 1.947 from ESR. 
Since merely scales with X: the underestimate is re- 
lated to the paramagnetic impurity which implies a re- 
duction in the contribution of the AgV0As04 phase. 
Assuming 7 % of the spin-i impurity, we arrive at 
g = 1.947 X a/1 - 0.07 ~ 1.88 in agreement with g ~ 
1.87 from the fit. We note that the impurity contribu- 
tion is not caused by a foreign crystalline phase (this 
is excluded by synchrotron XRD) and rather originates 
from AgV0As04 itself (e.g., from defects that break the 
spin chains). Exchange anisotropy of the Dzyaloshinsky- 
Moriya (DM) type can also be ruled out due to the crystal 
symmetry (see Sec. IVII) . 

The low-temperature magnetization curve of 
AgV0As04 (Fig. [3]) is characteristic of a spin-gap 
magnet. After subtracting the 7 % contribution of 
spin-i paramagnetic impurity approximated by a Bril- 
louin function,?3, we find that the magnetization (M) 
is close to zero below fioHc — 10 T, increases above 
10 T due to the closing of the spin gap, and saturates at 
IIqHs ~ 48.5 T. The bend at He directly measures the 
spin gap A = g^io^isHc/kB ^ 13 K3 

The comparison to the simulated curve for the alter- 
nating spin chain with J and J' from the susceptibility 
fit underscores two important differences. First, the fit- 
ted values of J and J' overestimate the spin gap and the 
saturation field. Second, the shape of the experimental 
curve is different from one expected for a purely ID sys- 
tem: the bend around 20 T is not observed, while the 
bend at 40 — 45 T, preceding the saturation, is notably 
smeared. Since we simulate the curve at a realistic tem- 
perature of T/ J = 0.035 (i.e., 1.4 K), the smearing of the 
bends cannot be ascribed to thermal fluctuations. The 
smooth shape of the experimental curve is rather caused 
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FIG. 3. (Color online) Top: high-field magnetization curve 
measured at 1.4 K and corrected for the contribution of para- 
magnetic impurities;^"^ solid and dashed lines show the fits of 
the ID (isolated alternating chain) and 2D (coupled alternat- 
ing chains) models, respectively. Experimental values of the 
magnetization (M) are scaled for the saturation magnetiza- 
tion (Ms). Bottom: the model of coupled alternating chains. 



by interchain couplings that modify the excitation spec- 
trum of the system. 

The 2D model of coupled alternating spin chains {J — 
J' — J±) improves the fit of the magnetization curve. 
We were able to impose either FM, AFM, or both FM 
and AFM interchain couplings to obtain nearly indistin- 
guishable fits (dashed line in Fig. |3] and Table ITV)) . The 
simulations are done in dimensionless units H/ (gfiBJ), 
hence both g and J values simply scale the field axis. 
Therefore, we fixed g = 1.947 according to ESR and ad- 
justed J. We find a higher J for FM interchain couplings, 
whereas AFM interchain couplings contribute to the sat- 
uration field and reduce J. However, none of the regimes 
fits to xiT) and M (H) simultaneously. For AFM inter- 
chain couplings, the fit of the high-field magnetization 
underestimates J compared to the susceptibility fit. For 
FM interchain couplings, the J values from the two fits 
nearly match, yet the fitted g-valne for x{T) is well below 
1.947 (or ~ 1.88 once the impurity contribution is taken 
into account, see above). The model with FM- and- AFM 
interchain couplings yields a reasonable value, although 
sizable discrepancies in the J value remain (see Table ITVT) . 

Presently, we are unable to find a model that will 
provide a consistent fit of the magnetic susceptibility 
and high-field magnetization. This situation is drasti- 
cally different from another alternating-chain-based sys- 
tem PbaVaOg (Ref. [Ml) where the introduction of FM 
interchain couplings readily resolves all the drawbacks of 
the ID model. The lack of the consistent fit likely re- 
sults from the complex and frustrated arrangement of 
the interchain couplings, as further confirmed by our mi- 
croscopic analysis in Sec. [Vj 



TABLE IV. Fitting of the magnetic susceptibility and high- 
field magnetization data with the model of coupled alternat- 
ing spin chains (bottom part of Fig. [3]): J_l = > (first 
line), J± — J'^ < (second line), and J± — —J'± (third 
line) . The magnetization curves were fitted with the constant 
g = 1.947 from ESR. 



Jx/J 


J'/J 


J (from x) 


g (from x) 


J (from 


M) 


0.15 


0.60 


42 


1.93 


34 




-0.27 


0.60 


41 


1.79 


40 




±0.15 


0.60 


41 


1.86 


37 






B. 


Electron 


spin resonance 






The 


ESR 


experiments 


on powder 


samples 


of 



AgV0As04 yield spectra with a shape typical for 
powder-averaged Lorentzian lines, see the inset of Fig. 2] 
The anisotropy of the g values corresponds to a reso- 
nance field anisotropy which is clearly resolved in the 
Q-band spectra. As shown by the solid fitting lines, 
the two components of a uniaxial powder fit are suf- 
ficient to describe the data (symbols) in a reasonable 
quality, despite the crystallographic symmetry is mon- 
oclinic. This fitting works well in the whole tempera- 
ture range under investigation. At the X-band frequency 
and T = 293.5 K, we obtained = 1.920 ± 0.002 and 
g± = 1.960 ± 0.002 which corresponds to an average g 
factor^ g = i(5|| + 2g±) = 1.947 ± 0.002. These values 
are similar to those found for V*"*" in other compounds 
with VOg octahedra, see, e.g., Ref. HI. The Q-band spec- 
tra can be described by the same (/-values within the ex- 
perimental error. Note that the difference between g± 
and g\\ implies the anisotropy of a single magnetic ion, 
and does not indicate exchange anisotropy. 

The temperature dependence of the integrated ESR 
absorption and x(T) (Fig. [2]) agree to each other (not 
shown), which confirms the validity of the temperature 
dependencies of the X-band ESR parameters shown in 
the main panels of Fig. 01 Both g- values remain constant 
in a wide temperature range down to 35 K and 22 K for 
(/|| and g±, respectively. The linewidth Ai3||^(T), de- 
scribing the half-width at half-maximum of the ESR ab- 
sorption, shows a characteristic behavior as observed in 
many other low-dimensional spin systemsi^"— It steeply 
increases below 50 K indicating the onset of spin corre- 
lations at low temperatures. Above 50 K, the anisotropy 
of the linewidth appears to be temperature-dependent: 
Ai3||(T) shows a shallow minimum at around 210 K 
whereas AB± continues to decrease with increasing tem- 
perature. The weak anomaly near 50 K is likely related 
to uncertainties in the fitting procedure. 
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FIG. 4. (Color online) Temperature dependence of X-band 
ESR parameters: g factor (top panel) and ESR linewidth 
AB (bottom panel). The inset shows the room-temperature 
spectra (symbols) at two microwave frequencies, X-band (9.4 
GHz) and Q-band (34 GIfz), together with powder-averaged 
Lorentzian line fits (solid lines) that provided the axial- 
symmetric parameters gn± and ABi\±. 



V. BAND STRUCTURE AND MICROSCOPIC 
MODEL 

To understand the relation between the ID 
alternating-spin-chain model and the 3D crystal 
structure of AgV0As04, we perform band structure 
calculations and evaluate individual exchange couplings. 
The LDA energy spectrum (Fig. [S]) reveals an electronic 
structure, typical for compounds. The low- 

lying valence bands are predominantly formed by O 2p 
states. The bands around —3 eV represent the filled Ad 
shell of Ag+i (similar to Ref. [H). The bands at the 
Fermi level have V 3d origin and overlap with empty 
states of Ag (5s) and As (4s) above 3 eV. A closer 
examination reveals the t2g and eg crystal-field levels 
of V 3d shell below and above 1 eV, respectively. The 
gapless energy spectrum results from the well-known 
underestimate of correlation effects in LDA. 

The formation of the short V-0 bond typically leads 
to a crystal-field split of the t2g levels. Then, the half- 
filled orbital has the xy symmetry and lies in the plane 




Energy (eV) 

FIG. 5. (Color online) LDA density of states for AgV0As04. 
The Fermi level is at zero energy. 



perpendicular to the short bondi^ '^^i'^^ The electronic 
structure of AgV0As04 reminds of this scenario, yet 
the t2g levels are not fully split. The bands between 
—0.2 and 0.2 eV indeed have a dxy origin but they are 
strongly hybridized with the higher-lying dyz and d^z 
bands (Fig. IB]). This can be related to the octahedral 
coordination of vanadium and to the relatively short 
V-0(3) bond (2.14 A) opposite to the vanadyl bond, sim- 
ilar to PbzVaOg (Ref.|3l|). 

To include correlation effects and evaluate exchange 
couplings, we construct a three-orbital tight-binding 
model (Fig. [SJ based on Wannier functions with proper 
orbital characters;^ The reduction to the three-orbital 
model is justified by the low contribution of the Cg 
states to the exchangei^ The orbital energies of 0.05 eV 
{xy), 0.50 eV (xz), and 0.60 eV (yz) are consistent with 
our expectations regarding the lowest energy of the dxy 
level. In the correlated Mott-insulating state (as obtained 
from LSDA-I-J7 calculations), the dxy orbital is half- filled, 
and the exchange couplings are derived from the Kugel- 
Khomskii model)^!^ 



J 



Ati 



off 



E 



([/off + A„)([/eff + A„ - Jeff)' 



(3) 



where txy and t^y 



are the hoppings between the xy 
states and from the xy (half-filled) to a (empty) states, 
UoS and Jeff are the effective on-site Coulomb repulsion 
and Hund's coupling in V 3c? bands, respectively, and 
Aa is the crystal-field splitting between the xy and a 
states. The first term in Eq. ^ is an AFM coupling due 
to hoppings between the half-filled orbitals. The second 
term stands for a FM coupling that arises from hoppings 
to empty orbitals followed by the Hund's coupling on 
the vanadium site. Note also that the system is in the 
strongly correlated regime {ti <C Ucs). 

To derive J's, we take UeS = 4 eV and Jeff — 1 eV, 
as established for other V"*"* compounds. ■^A"^'^ The re- 
sulting exchanges are listed in Table |Vl The leading 
couplings J and J' are AFM and run along [110] or 
[110]. This leads to alternating spin chains that align 
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TABLE V. Interatomic distances and exchange integrals 
(in K) calculated within the model approach {UeS = 4 eV, 
Jeff — 1 eV) and the supercell approach (LSDA+f/, AMF, 
I7d = 4 eV, Jd = 1 eV). 



Distance 


jAPM 


J™ J 


J 


(A) 


Model approach 


LSDA+17 


J 5.59 


88 


88 


47 


J' 5.56 


59 


-9 50 


47 


Jc 3.64 





-4 -4 


-9 


Ja 6.12 


11 


11 


6 


TABLE VI. Exchange integrals (in K) calculated via the su- 
percell approach with different values of LSDA+(7 Coulomb 
repulsion parameter Ud (in eV) and different DCC schemes, 
see text for details. 


Ud J 


J' 


Jc 


Ja 


3 62 


61 


-2 


7 AMF 


4 47 


47 


-9 


6 AMF 


5 35 


38 


-24 


5 AMF 


4 41 


36 


-8 


4 FLL 



with different crystallographic directions to form an un- 
usual crossing pattern (middle panel of Fig.[T]). The cou- 
plings between the spin chains, FM Jc (along the struc- 
tural chains) and AFM Ja (approximately along [102]), 
complete our model to a peculiar 3D spin lattice. Further 
couplings are below 1 K and can be neglected. 

LSDA-f J7 supercell calculations (last column of Ta- 
ble |V| confirm the results of the model analysis. We find 
leading AFM couplings J and J' along with weaker in- 
terchain couplings, FM Jc and AFM Ja- Note, however, 
that the alternation ratio can not be evaluated precisely. 
The model approach yields a = J' / J = 0.57 in remark- 
able agreement with the experimental a = 0.60 — 0.65, 
yet the absolute values of J and J' are overestimated by 
a factor of two (e.g., J = 88 K in DFT vs. J ~ 40 K in 
the experiment). LSDA+U better evaluates the absolute 
values but overestimates the a ratio. 

To get additional insight regarding the influence of 
computational parameters on individual exchange cou- 
plings, we calculated exchange integrals for different val- 
ues of Ud, the on-site Coulomb repulsion parameter of 
the LSDA+C/ method (Table EH). We also applied differ- 
ent double-counting correction (DCC) schemes, around- 
mean- field (AMF) and fully- localized-limit (FLL), since 
the choice between the two is known to affect exchange 
couplings in Cu"*"^ compounds.— "^"^ Surprisingly, the ex- 
change couplings in AgV0As04 are largely insensitive to 
the DCC, whereas the change in Ud leads to a typical 
evolution of J's (see, e.g., Ref. The a ~ 1 ratio 

is obtained for a range oi Ud values, hence this result 
should be taken as an intrinsic drawback of LSDA+C/ 
calculations for AgV0As04. Overall, we find that both 




rZ T R A Y M X V 



FIG. 6. (Color onhne) LDA band structure of AgV0As04 
(thin light lines) and the fit with the tight-binding model 
(thick dark lines). The fc-space path is defined as fol- 
lows: r(0,0,0), Z(0, 0,0.5), r(0.5, 0,0.5), i?(0.5,0.5,0.5), 
A(0, 0.5, 0.5), y(0, 0.5, 0), M(0.5, 0.5, 0), X(0.5, 0, 0) with the 
coordinates given in units of the respective reciprocal lattice 
parameters. The Fermi level is at zero energy. 

the model and supercell approaches experience certain 
difficulties in the precise evaluation of the exchange cou- 
plings. The combination of different methods allows to 
derive the complete spin model, and facilitates the micro- 
scopic understanding of the material, whereas the exten- 
sive analysis of the experimental data (Sec. IIVP justifies 
our assignment of AgV0As04 to the model of coupled 
alternating spin chains. 



VI. DISCUSSION AND SUMMARY 

The experimental results for the magnetic suscepti- 
bility and high-field magnetization of AgV0As04 can 
be qualitatively understood within the bond-alternating 
spin chain model. The alternation of the exchange cou- 
plings J and J' opens the spin gap A ~ 13 K, and leads 
to the likely singlet ground state (no long-range order in 
zero field) . The microscopic origin of J and J' is rational- 
ized by DFT. Inversion centers of the P2i/c space group 
relate the neighboring vanadium atoms within the spin 
chains, thereby DM couplings on the J and J' bonds are 
forbidden by symmetry. Despite the weak anisotropy of 
the g-tensor, the exchange couplings in AgV0As04 are, 
to a good approximation, isotropic. 

The alternating-chain magnetic behavior is rather 
common for low-dimensional magnets. It has been previ- 
ously observed in PbsVaOg (Refs.EUndlil), Cu(N03)2- 
2.5H2O (Ref. 0), CUWO4 (RefHli), and other com- 
pounds. A distinctive feature of AgV0As04 is the un- 
usual regime of the interchain couplings. A close ex- 
amination of the spin lattice shows that the interchain 
couplings Ja and Jc are frustrated. The four-bond unit 
depicted by the thick line in Fig. [T] comprises one FM 
and three AFM couplings, which can not be satisfied si- 
multaneously (right panel of Fig. [T|). 

Magnetic frustration has a dramatic effect on thermo- 
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dynamic properties and ground state of spin-chain mag- 
nets. For example, frustrated interchain couphngs in- 
duce an incommensurate magnetic structure in CS2CUCI4 
(Ref. "46^ and a ^-magnetization plateau in Cs2CuBr4 
(Ref. 47). Ising-like systems, such as Ca3Co206 and 
CaaCoRhOg, show an even more intricate behavior re- 
lated to the formation of a partially ordered AFM 
stateJ^!^ The frustrated nature of interchain cou- 
plings could similarly distinguish AgV0As04 from other 
alternating-spin-chain magnets. The leading energy scale 
is, however, determined by the alternating J and J' 
bonds, and dictates the spin gap formation. Above the 
first critical field iJc, external magnetic field closes the 
spin gap, thereby a long-range magnetic order (BEC 
of magnons) is possiblciii^i^ Since the frustration im- 
pedes the long-range order, one might expect the par- 
tial suppression of the BEC. The accessible critical fields 
iHoHc ^ 10 T, unHs ~ 48.5 T) make AgV0As04 an in- 
teresting compound for detailed high-field studies apply- 
ing specific heat, magnetocaloric effect, ESR, and nuclear 
magnetic resonance measurements. 

Turning to the structural aspects of the proposed 
spin model, we note that the frustration scenario in 
AgV0As04 is very similar to one in M(VO)2(P04)2 with 
M — Ca and SiM- In both cases, most of the couplings are 
long-range and AFM, yet the presence of the short-range 
FM coupling (Jc in AgV0As04) between the corner- 
sharing VOg octahedra makes the system frustrated. The 
magnitude of the frustration is controlled by the ratios 
of J's. In M(VO)2(P04)2, four comparable couplings in 
the frustrated unit lead to the strong frustration on a 
3D spin lattice, and the enhancement of quantum fluctu- 
ations manifests itself as the suppression of the specific 
heat maximum. In AgV0As04, J and J' dominate over 
Ja and Jc, hence the magnetic behavior is quasi-lD, with 
the spin gap imposed by the alternating spin chains (sim- 
ilar to non- frustrated systems). The effects of the frus- 
tration should be primarily expected in high magnetic 
fields that close the spin gap and induce the long-range 
order. 

The peculiar superexchange regime of AgV0As04 is 
readily understood from the examination of Wannier 
functions (Fig. [T]) . Owing to the symmetry of the half- 
filled vanadium orbital, the Wannier functions do not 
comprise the orbitals of the 0(3) atom, and render the 
V-O-V superexchange pathway (Jc) inactive. Indeed, 
we find a weak FM exchange for Jc, despite the shortest 
V-V separation at this pathway. The leading exchanges 
J and J' run via longer pathways with double bridges of 
ASO4 tetrahedra, since the edges of the tetrahedra pro- 
vide a suitable overlap of the oxygen orbitals. Finally, 
the weaker AFM coupling Ja can also be traced to a fa- 
vorable overlap of the oxygen orbitals on the edge of the 
tetrahedron (a single bridge). All these findings conform 
to previous, both empirical^ and microscopio^s^ con- 
siderations regarding the superexchange in systems with 
tetrahedral polyanions: (i) the possible superexchange 
pathways are determined by the symmetry of the mag- 




FIG. 7. (Color online) Wannier functions for the superex- 
change pathways Jc (left), J and J' (top right), and Ja (bot- 
tom right). Thick (orange) lines denote short V-0(3) bonds. 



netic orbital; (ii) the tetrahedral anionic groups lead to 
short 0-0 distances and induce suitable overlaps of oxy- 
gen orbitals; (iii) double anionic bridges are more effi- 
cient than single bridges, owing to a different geometry 
and a larger number of 0-0 contacts. We expect that 
these simple arguments will facilitate further search for 
unusual spin lattices in V*"^ phosphates, arsenates, and 
other compounds with tetrahedral anions. 

In summary, we have studied the crystal structure and 
magnetic properties of AgV0As04, and proposed a mi- 
croscopic description of this compound. This system ba- 
sically conforms to the alternating spin chain model, al- 
though the quantitative understanding of the magnetic 
behavior requires an extension of the model toward a 3D 
spin lattice with frustrated couplings between alternat- 
ing spin chains. Our magnetic susceptibility and high- 
field magnetization measurements establish a spin gap of 
about 13 K and a saturation field of 48.5 T. These es- 
timates suggest the feasibility of further high-field mea- 
surements that should explore the influence of frustra- 
tion on the behavior of dimer-based quantum magnets. 
The study of the Bose-Einstein condensation phenom- 
ena in AgV0As04 might be particularly enlightening and 
calls for extensive experimental investigation of this com- 
pound. 
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FIG. SI. Crystal structure refinement at 300 K: experimental (symbols), calculated (yellow line), and difference (black line) 
patterns. Ticks show the reflections of AgV0As04 (lower line) and the Ag impurity (upper line). 



1.0 
0.8 



^0.6 



g 0.4 
e 

0.2 
0.0 

0.0 



0.2- 

0.1 

0.0- 
0.0- 



M — ^- 



9 10 11 

JijiliuJi,ilftJi,JijJj'bjulJlLjl*L 



12 



13 



14 



1 III! lllllllllllll 

— ,—uJ\ 


II II 




II II 













10 15 20 

Scattering angle 20 (deg) 



25 



30 



FIG. S2. Same as Fig. ISll for the structure refinement at 20 K. 
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Temperature (K) 



FIG. S3. Magnetic susceptibility of AgV0As04 measured in the applied fields of 0.5, 1, 3, and 5 T. The arrow shows the 
reduction in the low-temperature upturn upon the increase in the magnetic field H . 




Field (T) 

FIG. S4. Experimental high-field magnetization curve measured at 1.4 K (circles), the contribution of the paramagnetic 
impurity (dashed line), and the intrinsic magnetization (solid line) of AgV0As04 obtained by the subtraction. The inset shows 
the field derivative of the experimental magnetization curve along with the critical fields He (closing of the spin gap) and Ha 
(saturation) . 



